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The evolution of tetragonality with thickness has been probed in epitaxial c-axis oriented PbTiOa 
films with thicknesses ranging from 500 down to 24 A. High resolution x-ray pointed out a systematic 
decrease of the c-ax;is lattice parameter with decreasing film thickness below 200 A. Using a first- 
principles model Hamiltonian approach, the decrease in tetragonality is related to a reduction of 
the polarization attributed to the presence of a residual unscreened depolarizing field. It is shown 
that films below 50 A display a significantly reduced polarization but still remain ferroelectric. 
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Since its discovery in 1920 by Valasek ferroelec- 
tricity has attracted considerable interest from a fun- 
damental point of view and because of its wide range 
of potential applications. Ferroelectricity has been his- 
torically seen as a collective phenomenon [3| requiring a 
relatively large critical volume of aligned dipoles. Calcu- 
lations within the Devonshire-Ginzburg-Landau (DGL) 
theory had indicated that, as the physical dimensions of 
a ferroelectric material are reduced, the stability of the 
ferroelectric state is altered, leading to a relatively large 
critical size below which ferroelectricity is suppressed in 
small particles and thin films 0, Q- For instance, in the 
case of Pb(Zro. 5X10.5)03 thin films, a critical thickness 
of ~200 A had been predicted at room temperature 0. 
These predictions appeared to agree with then-current 
experiments Recent results, however, suggest a much 
smaller critical size, with ferroelectricity detected in poly- 
mer films down to 10 A j™ and in perovskite films down 
to 40 A (10 unit cells) X-ray synchrotron studies 
have revealed periodic 180° stripe domains in 12 to 420 
A thick epitaxial films of PbTiOa (PTO) grown on insu- 
lating SrTiOs (STO) substrates [101. On the theoretical 
side, atomistic simulations have emphasized the predom- 
inant role of electrostatic boundary conditions in deter- 
minin g fe rroelectricity in very thin films. Ghosez and 
Rabe [l3|, and Meyer and Vanderbilt jldj showed that, 
under perfect screening of the depolarizing field, ultra- 
thin stress- free perovskite slabs can sustain a polarization 
perpendicular to the surface. However, perfect screening 
is not achieved in usual ferroelectric capacitors. Batra 
et al. [T^ showed that, under short-circuit boundary 
conditions, the incomplete screening of the depolarizing 
field resulting from the finite screening length of the metal 
can substantially affect the ferroelectric properties [sof . 
Consequences of imperfect screening on the coercive field 



* Electronic address: [Celine. Lichtensteiger@physics.unige. 

t Present address : Department of Physics and Astronomy, Rutgers 

University, Piscataway, NJ 08854-8019, USA 



have recently been studied in detail [ig. Also, first prin- 
ciples calculations allowed the amplitude of the depo- 
larizing field for BaTiOa ultrathin films between metallic 
SrRuOa electrodes in short-circuit conditions to be quan- 
tified and a critical thickness of ~24 A (6 unit cells) in 
such structures to be predicted [T^ . 

These recent experimental and theoretical results, at 
odds with the former well established belief, clearly show 
that additional studies are crucial to demonstrate exper- 
imentally the key role of screening and confirm ferroelec- 
tricity at the nanoscale ,18.] . In this Letter we report 
on combined experimental and theoretical investigations 
of tetragonality in a series of epitaxial c-axis oriented 
films of PTO grown onto metallic Nb-doped STO sub- 
strates. X-ray analyses show that the tetragonality pro- 
gressively decreases below 200 A. Using a first-principles- 
based model Hamiltonian approach, we relate this low- 
ering of tetragonality to a lowering of the spontaneous 
polarization of the films, due to a residual unscreened 
depolarizing field. Our results show that ultrathin PTO 
films on Nb-STO display a significantly reduced sponta- 
neous polarization but still remain ferroelectric below 50 

A. 

Two series of PTO films were grown onto metallic (001) 
0.5 wt% Nb-STO substrates by off-axis radio-frequency 
magnetron sputtering |0, |23| , with 60W and 40W ap- 
plied to the Pbi.iTiOs target, corresponding to growth 
rates of 280 A/h and 110 A/h respectively [3ll |. 

Room temperature x-ray measurements, using a 
Philips X'Pert High Resolution diffractometer, allowed 
us to determine precisely the epitaxy of the films, their 
thickness and their c-axis parameter. The 9-29 diffrac- 
tograms, revealing only (00^) reflections, demonstrate 
that the films are purely c-axis with the polarization 
normal to the film surface. Pole-figures confirmed the 
expected epitaxial "cube on cube" growth and that the 
films were tetragonal. High angle finite size oscillations 
and low angle reflcctometry (Fig. top) allowed us to 
determine, through simulations, the number of planes in- 
volved in the diffraction and thus the film thickness 
as well as the deposition rate, even for films down to 24 
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FIG. 1: Top: 61-261 difltractogram around the (001) diflraction 
peak for a 396 A sample (thick line) and simulation (thin line) . 
Inset: low angle 6'-26' diffractogram for the same sample. Mid- 
dle: 9-26 diffractogram revealing (GO?) reflections from 1=1 to 
5 for the same 396 A thick fihn. Bottom: 6-26 diffractograms 
around the (002) diffraction peak for three films of different 
thicknesses: 396, 132 and 73 A respectively, allowing the de- 
crease of the c-axis with thickness to be clearly seen. 



A. To precisely determine the c-axis parameter, we per- 
formed x-ray diffraction for the (OOZ), l—l to 5 reflections 
as shown in Fig. ^ (middle) . Since at room temperature 
the a-axis latti ce p arameter of ferroelectric bulk PTO is 
3.902-3.904 A HI El, very close to the 3.905 A of the 
STO, the films are expected to be coherent, with their 
a-axis equal to the STO lattice parameter. Grazing inci- 
dence diffraction on a ~86 A thin film confirmed this pic- 
ture, displaying a unique (200) reflection. Measurements 
of the PTO (101) reflection for the thickest films (-504 
A and ^^396 A) gave an estimation of a = 3.90 ± 0.01 A. 
For these films, in-plane strain relaxation would lead to 
a maximum change in the c-axis length of 0.006 A, which 
is ~15 times smaller than the changes discussed below. 
The a-axis value used later to calculate the c/a ratio is 
thus taken as constant and equal to 3.905 A. 

To probe finite size effects in thin films, a key require- 
ment is to have materials with smooth surfaces, and 
therefore a well-defined thickness. Atomic force micro- 
scope (AFM) topographic measurements were performed 
on all the samples, showing that the films are essen- 
tially atomically smooth with a root-mean-square (RMS) 
roughness between 2 and 6 A over 10x10 ^m? areas. 
Fig. [3 (left) shows a representative topographic image 
obtained on a 33 A thick film. The vertical scale used 




FIG. 2: Left: AFM topography for a 33 A thick film. The 
RMS roughness is ~ 3 A over the scanned area. The verti- 
cal scale is the film thickness. Right: Piezoresponse signal 
obtained after writing ten stripes using alternate -I-12V and 
-12V voltages applied to the tip. 



in this 3D representation is equal to the thickness of the 
sample, allowing a comparison of the roughness of the 
surface with the total film thickness. As can be seen from 
the data, the average thickness remains well defined, even 
for ultrathin films. 

Finally, we used the piezoresponse mode of the AFM 
to probe the domain structure of different films |3j IH ■ 
Ten stripes where drawn using alternate -I- 12V and -12V 
voltages applied to the tip over a 1.6x1.6 ^w? area. The 
3.2x1.6 /im^ background piezoresponse signal was then 
compared to the signal from the stripes (Fig. [21 (right)) 
and found to be equal to the signal obtained for the +12 V 
written stripes suggesting a single polarization in the as 
grown sample. All the investigated films were found to 
be "monodomain-like" over the studied areas, typically 
a few /xm^. 

We now return to the relation between the c-axis pa- 
rameter and the film thickness. Fig. (bottom) is a 
blow up of diffractograms around the (002) diffraction 
peaks and allows the change in lattice parameter with 
film thickness to be directly seen for three samples. Fig.|3| 
shows the central result of the paper : the film tetrago- 
nality, i.e. the c/a ratio, is plotted as a function of film 
thickness for the two series of samples (top). As can be 
seen, the data for both series collapse and the c/a ratio 
decreases very substantially for films thinner than 200 A. 

We note that a similar reduction of tetragonality has 
recently been observed by Tybell [2^. It has also been 
checked experimentally that the c-axis values did not 
change after the deposition of a gold electrode and short- 
ening of the gold electrode and the metallic substrate 
(insuring short circuit conditions). 

The observed reduction of c/a shown on Fig. O (top) 
cannot be attributed to a change in the a-axis lattice pa- 
rameter since it is independent of the film thickness. In- 
stead, the decrease of c must be related to a concomitant 
reduction of the spontaneous polarization through the 
polarization-strain coupling that is known to be partic- 
ularly large in PTO j23|. The evolution of c/a therefore 
is a signature of the progressive suppression of ferroelec- 
tricity in ultrathin films. 

This can be highlighted theoretically by extending the 
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FIG. 3: Evolution of the c/a ratio with the film thickness. 
Top: experimental results for the 1st series (circles), the 2nd 
series (squares), and the 1st series with a gold top electrode 
(crosses) ; the dashed line is the phenomenological theory pre- 
diction (see main text) supposing a ratio between the extrapo- 
lation and the correlation length S/^ = 1.41 [2^ : the solid line 
is the model Hamiltonian prediction for Ae//=0.12 A, rescaled 
to give a maximum tetragonality in agreement with the ex- 
perimental data. Bottom: results from the model Ifamilto- 
nian calculations for Ae//=0.23 A (dashed line), Ae//=0.12 A 
(solid line) and Ae//=0.06 A (dotted line). Inset: thickness 
dependence of the spontaneous polarization V calculated from 
the model Hamiltonian for Ae//=0.12 A. 



first-principles effective Hamiltonian approach developed 
by Waghmare and Rabe for bulk PTO |23 to thin films. 
Within this approach, the energy is written as a low- 
order Taylor expansion of the bulk energy around the 
cubic phase within the restricted subspace spanned by 
(i) the ionic degree of freedom ^ associated to the lat- 
tice Wannier function of the soft phonon branch and (ii) 
the macroscopic strain e. The bulk Hamiltonian includes 
three terms : the double-well soft-mode energy, the elas- 
tic energy and a coupling term between ^ and e. All the 
bulk parameters are obtained by fitting first-principles 
results. 

InRef.lIl the suppression of the ferroelectricity in ul- 
trathin films between metallic electrodes in short-circuit 
was related to the incomplete screening of the depolar- 
izing field by the electrodes. It was also shown that the 
energetics of ultrathin films can be reasonably described 
using a simple model that corrects the bulk internal en- 
ergy at first order by a term taking into account the en- 
ergy associated to the coupling between the polarization 
V and the unscreened depolarizing field Sd- Transposing 
this to the effective Hamiltonian approach allows us to 
write : 

nT[C,e]=H^J^[te]-£d-V (1) 

where 7^^^"^ and Tioff" ^^e the model Hamiltonians for 
bulk and thin film respectively. For TC^g^^, we keep the 



form and parameters fitted at the experimental volume 
as reported in Ref. Therefore, the only finite-size 
correction included in arises from the electrostatic 

energy related to £d- 

Under short circuit conditions, £d is related to the po- 
tential drop Ay at each metal/ferroelectric interface. As- 
suming two similar interfaces with the top and bottom 
electrodes, then £d — —2AV/d, where d is the thick- 
ness of the film. The potential drop was attributed to 
finite dipole densities at the interfaces and was shown 
to evolve linearly with the spontaneous polarization : 
AV — (Ae///eo) -V. The parameter Ag// has the dimen- 
sion of a length and will be referred to as the effective 
screening length of the system. 

Supposing monodomain films polarized along z, in 
agreement with the piezoresponse measurements, the 
macroscopic polarization is homogeneous and directly 
linked to i {Vz = Z*^^/^.^ ^ where Z* is the soft 
mode effective charge, £^z is the z-componcnt of ^, and 
ilo is the unit cell volume) so that we finally obtain : 

-£d-Vz^ {2X,ffZ*ynleod)ez (2) 

This energy scales with X^ff/d and is positive meaning 
that the effect of the depolarizing field is to suppress the 
ferroelectric instability through a renormalization of the 
quadratic term of the effective Hamiltonian. 

The previous model is now applied to thin films of PTO 
on Nb-STO. Perfect pseudomorphic thin films on top of 
a cubic substrate are considered and the in-plane strains 
(exx = eyy = {asTO - apTo)/apTO with asTO = 3.905 
A and apTO = 3.969 A ||3^) are fixed throughout the 
structure independently of the film thickness, in order 
to constrain the a-axis lattice constant imposed by the 
STO substrate. The energy (Eq.^ is then minimized for 
different thicknesses in terms of (supposed uniform and 
perpendicular to the film) and e^z ■ From the values of 
and Czz we deduce V and c/a. All the simulations have 
been carried out at T = 0. 

First-principles results for the SrRuOa/BaTiOs inter- 
face yield Ag// = 0.23 A. Because the screening might 
be slightly different for the present system composed of 
distinct interfaces, the theoretical results are reported 
in Fig. 121 (bottom) for slightly different values of Xefj- 
The model predicts a critical thickness below which the 
spontaneous polarization vanishes and the c/a ratio of 
the resulting paraelectric phase saturates at 1.03, as a 
result of the mechanical constraint imposed by the sub- 
strate. Above the critical thickness, the spontaneous po- 
larization gradually increases up to the bulk value as 
does the c-axis lattice parameter, as a consequence of 
the polarization-strain coupling. 

The model Hamiltonian of Waghmare and Rabe, al- 
though appropriately describing PTO, is known to over- 
estimate the polarization-strain coupling. At the bulk 
level, the model predicts c/a= 1.09 while the experimen- 
tal value is equal to '^1.06 [l^. In order to get rid of this 
bulk overestimation, and for a direct comparison of the 
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theoretical and experimental evolution of c/a, the theo- 
retical curves have been renormalized to give a tetrago- 
nality of 1.068 at 500 A, in agreement with the experi- 
mental data. Only the strength of the polarization-strain 
coupling must be rescaled while the c/a value of the para- 
electric phase, a priori properly predicted through the 
elastic constants, can be kept unchanged. 

Looking at the experimental points on Fig. |3| (top), 
both the range of thicknesses at which the c/a ratio starts 
to decrease and the shape of the evolution agree with 
the prediction of the model Hamiltonian calculations for 
Ae//=0.12 A. This supports an incomplete screening of 
the depolarizing field as the driving force for a global re- 
duction of the polarization in perovskite ultrathin films. 
In contrast, the DGL phenomenological theory including 
only an intrinsic suppression of ferroelectricity at the sur- 
face through the so-called extrapolation length parame- 
ter S .33.] predicts a much sharper decay with no substan- 
tial decrease of polarization predicted above 50 A 27] . 

The value of Ae// (0.12 A) used in Fig. 01 (top) is 
smaller than in the case of the BaTiOs/SrRuOa inter- 
face ^3 (0.23 A) (albeit in the same order of magni- 
tude). This might suggest a better screening for the 
present system but might also be partly attributed to 
an overestimate of the theoretical critical thickness due 
to the simplicity of the model Hamiltonian "si"] and to 
the fact that the simulations were performed at T=0. 

Importantly, the thinnest films have a much higher 
tetragonality than the value of 1.03 expected for the para- 



electric phase from the macroscopic theory of elasticity. 
No saturation of c/a, the signature of a complete sup- 
pression of ferroelectricity, was observed, clearly implying 
that films much thinner than 50 A are still ferroelectric. 

In conclusion, we found that the c-axis parameter of 
PTO films decreases substantially below 200 A. A first- 
principles effective Hamiltonian approach allowed us to 
establish that the lowering of c/a is related to a progres- 
sive reduction of the polarization due to an imperfect 
screening of the depolarizing field. Although their polar- 
ization is significantly reduced, the fact that no satura- 
tion of the c-axis to its paraelectric value was found down 
to very thin films (24 A) demonstrates that PTO films 
below 50 A remain ferroelectric at room temperature. 
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